Quantum phase transition in the one-dimensional extended Peierls-Hubbard model 



H. Benthien 1 , F. H. L. Essler 2 and A. Grage 1 
1 Fachbereich Physik, Philipps-Universitiit Marburg, D- 3 5032 Marburg, Germany 
2 The Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford 0X1 3NP, UK 

We consider the one-dimensional extended Hubbard model in the presence of an explicit dimerization S. 
For a sufficiently strong nearest neighbour repulsion we establish the existence of a quantum phase transition 
between a mixed bond-order wave and charge-density wave phase from a pure bond-order wave phase. This 
phase transition is in the universality class of the two-dimensional Ising model. 
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I. INTRODUCTION 

It is well known that quasi one-dimensional electron sys- 
tems exhibit a "Peierls instability" towards the formation of 
a dimerized insulating ground state |]|. In the absence of 
electron-electron interactions the low temperature phase of 
such systems is described in terms of a Peierls insulator with 
gapped electron and hole quasiparticle excitations. On the 
other hand, it is well known that undimerized interacting one- 
dimensional electron systems are either correlated (Mott) in- 
sulators or Luttinger liquids, see e.g. [2l|3j]. A characteristic 
feature of these states is that low-lying excitations are collec- 
tive modes of charge and spin degrees of freedom respectively. 
An interesting question is then what happens in the case when 
strong electron-electron interactions compete with the Peierls 
distortion. A simple model having all the necessary ingredi- 
ents to study this question is the one-dimensional extended 
Hubbard model with an explicit dimerization S. Its Hamilto- 
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where c[ a creates an electron with spin a =f, \ in a Wannier 
orbital centered around site I, and we have ni t(T — c~l a ci, a , 
ni = hi^ + hi t i, U is the on-site and V the next-neighbour 
Coulomb interaction. Since we are interested in the half-filled 
case only, the number of electrons N equals the number of 
lattice sites L. In what follows we will consider the charge- 
density wave (CDW) and bond-order wave (BOW) order pa- 
rameters 

meow = -J2(-l) l (ni-l) , (2) 
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m BO w = -^(-1)' (c^Q+i^+h.c.) . (3) 
The model Q has previously been studied in various lim- 



iting cases. The infinite-?/ limit was studied in |4|]. The low- 
lying excitations in this limit are chargeless spin triplet and 
spin singlet excitations, which can be understood in terms of 
a spin-Peierls Hamiltonian, see e.g. S 0- The effects 
of weak electron-electron interactions in a Peierls insulator 
were investigated in [6] by perturbative methods. Most impor- 
tantly for our purposes, the weak-coupling regime U, \ V\ < t 
was studied in Refs Qlal- From the structure of the classi- 
cal ground state of the bosonized low-energy effective Hamil- 
tonian Tsuchiizu and Furusaki showed that increasing the 
dimerization 5 from zero for a sufficiently large V drives the 
system through a quantum critical point that was argued to be 
in the universality class of the two-dimensional Ising model 
121. The mechanism underlying this transition is very similar 
to the one exhibited in iflOt ll ill The purpose of the present 
work is to verify the prediction 1 7] of an Ising critical point by 
means of numerical methods. The outline of the paper is as 
follows. In section||l]we consider the field theory limit of the 
model ([0 and review and extend the results of Ref. 0] that 
suggest the existence of an Ising transition. In section Hm we 
use numerical techniques to establish that there is indeed an 
Ising transition in the lattice model Q. 



II. FIELD THEORY LIMIT 

The weak-coupling limit U, V <C t of the model Q is 
amenable to a field theory analysis. The low-energy regime 
can be described by linearising the non-interacting Fermi 
spectrum around the Fermi -points ±7r/2ao, where ao is the 
lattice spacing. Applying a bosonization scheme 0, fp3l in 
terms of two Bose fields $ c and $ s corresponding to collec- 
tive charge and spin degrees of freedom respectively we arrive 
at the following form of the low-energy Hamiltonian | TlM f!^ 
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Here wp = 2tao is the Fermi velocity of the noninteracting 
theory and 8 C S are the fields dual to <J? CiS . The bare values of 
the couplings /3 C , g CtS and charge and spin velocities v c . s are 
related to the parameters of the lattice model as follows (see 
e.g. M) 



/3c = 




In H s we have kept the quadratic term in derivatives in the 
interaction part in order to emphasize the SU(2) symmetry of 
TC S - It is of course possible to absorb this term into the kinetic 
piece of the Hamiltonian through a rescaling of <£> s — » /3 S $ S , 
S — > /3^ 1 s . The effective low-energy model consists of 
two coupled sine-Gordon models (sGM) and cannot be solved 
exactly. The bosonized expressions for the order parameters 
are 
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From now on we restrict our analysis to the regime 2V > U, 
which in the absence of dimerization 8 = corresponds to the 
charge-density wave (CDW) regime. In this regime (for 5 = 
0) the Hamiltonian reduces to the well-known description of 
the extended Hubbard model in terms of two sGMs llal . The 
charge sector is described by a sGM with coupling constant 
(3 C < 1. As a result the charge sector is gapped and the gap 
scales like 
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Excitations in the charge sector are scattering states of gapped, 
spinless "(anti)holons" carrying charge =pe. The interaction 
of spin currents in the spin sector is marginally relevant and 
opens up a spectral gap (see e.g. 0E3l), which scales like 
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The elementary excitations are charge neutral spin-i spinon 
excitations with a spectral gap given by 1141 . 



A. Classical Ground State 



The qualitative behaviour of the field theory (@J can be de- 
termined by considering the classical limit f7l 1 1 Oil . The effec- 
tive potential is given by 
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where $ c = /3 C $ C , $ s = /3 s <f> s , In the charge-density wave 
regime U < 2V we have g c < 0, g s < and S > 0. The 
structure of the local minima of U c s then depends on the value 
of 5 as follows. 

1 . In the pure charge-density wave phase (5 = the min- 
ima are at {<J> C , ^s} = {(2n + l)7r, 2kir}, where n and k are 
integers. The pinning of the fields at these values implies the 
presence of a nonzero CDW order parameter (jticdw) 0. 

2. In the interval < 5 < 5* = \g c \/(2n) the num- 
ber of minima is unchanged, but their positions are shifted 
to {27r(2n + 1) ± (f> , 27r(2fc)}, {2?r(2n) ± <j) , 2it(2k + 1)} 
where <fio — 2 arccos(27r^/|5 c |). There are two nonzero order 
parameters 
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For small dimerization i « 1 the CDW order parameter is 
large compared to the BOW order parameter and decreases 

quadratically (meow) (pr) ■ 

3. At a critical dimerization 6 = 6* the adjacent minima 
that were moving towards each other merge. Like in the dou- 
ble sine-Gordon case 1 16] the analogy with the (^-description 
of the Ising model suggests that a quantum phase transition in 
the Ising universality class takes place at S*. 

4. For large dimerizations 6 > 6* the positions of the 
minima of U c g are independent of 6 and occur at {2-7r(2n + 
l),27r(2fc)}, {27r(2n),27r(2fc + 1)}. These minima are lo- 
cated at the same positions as in the pure Peierls insulator. As 
a result the CDW order parameter now vanishes, whereas the 
BOW order parameter stays finite 



(tocdw) = 0, (m BO w) ^ 0. 



B. Perturbation Theory in 8 
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a. Form Factor Perturbation Theory In the absence of 
dimerization 6 = the field theory (|5}, is integrable. Us- 
ing the knowledge of matrix elements (form factors) of oper- 
ators in this integrable theory JTtL [Tel . the effects of 8 > 
can be studied by form factor perturbation theory dlH. In 
the CDW regime 2V > U we are dealing with a fully massive 
quantum field theory. The changes in the holon and spinon 
gaps to first order in 5 are [ 19] 
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Here ZJ (9) is a Faddeev-Zamolodchikov operator creating a 
spinon with momentum A s sin h{9)/v s . Similarly Z\(B) cre- 
ates a holon with momentum A c sinh(0) / v c . The form factors 
in Jl 8i have been calculated in [17] and substituting them into 
(EH) we find 
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where £ = 0^/(1 — 01) and S r is the renormalized coupling 
constant at a scale set by the gap for 6 = 0. The fact that the 
corrections to the gaps are finite is an indication that perturba- 
tion theory in 8 is well defined at least for very small values 
of 5. 

b. Order Parameters for small 5 Both the classical anal- 
ysis and the formfactor perturbation theory considerations 
suggest that for sufficiently small 5 the Peierls term can be 
treated in perturbation theory around the gapped CDW phase. 
The same then ought to be the case for the lattice model itself. 
Let us then work with the lattice Hamiltonian and represent 
it as H = Hq + 6 Hi. The unperturbed extended Hubbard 
Hamiltonian H Q is invariant under the following discrete sym- 
metry transformation 11211 



(20) 



which is a combination of a particle-hole transformation and 
a translation by one site. It is straightforward to see that Hi is 
odd under this transformation 



UHiU^ = -Hi 



(21) 



On the other hand, the order parameters are odd and even re- 
spectively 



Um CTSW U^ = m CDW . 



(22) 



It follows from (12 U and (1221 that the perturbative expansions 
of the order parameters in powers of S are of the form 



(tocdw) = a n S 2n , 
n=0 



(23) 
(24) 



n=0 



The perturbative results for the order parameters, which obvi- 
ously can be derived in the field theory limit as well, show the 
same kind of dependence on S as the one obtained from the 
analysis of the classical ground state. 



C. Quantum Critical Point 

The 1-loop renormalization group analysis carried out in 
Ref.|7] suggests that in the regime we are interested in the 
spin degrees of freedom have a "large" gap and the low-energy 
effective Hamiltonian only involves the charge sector and is 
given by a two frequency sGM 
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where g*, g* s and j3* are renormalized c oup ling constants. It 
then follows from the analysis of II (All 1 111 611 that at some crit- 
ical value of 8 the charge sector undergoes a quantum phase 



transition in the universality class of the two-dimensional 
Ising model. At the transition point the charge degrees of 
freedom become gapless, while the spin degrees of freedom 
remain gapped. In the following section we carry out a nu- 
merical analysis of the underlying lattice model in order to 
assess the validity of this scenario. 



III. NUMERICAL RESULTS 

In this section we present numerical results for the CDW- 
and BOW order parameters, excitation gaps and critical expo- 
nents of the Hamiltonian Q obtained with the Density-Matrix 
Renormalisation Group (DMRG) O. DMRG is known to 
give excellent results for ground state expectation values and 
energies of one-dimensional lattice Hamiltonians and has be- 
come a standard method in the field. We show that the quali- 
tative results of the preceding field theoretical analysis are re- 
produced with DMRG and make a strong case that the quan- 
tum phase transition belongs to the Ising universality class. 
We note that the parameters U, V chosen in the numerical 
analysis are such that the field theory description @ is no 
longer quantitatively valid as can be seen from the fact that the 
numerically determined gaps are no longer small compared to 
the bandwidth it, which serves as the cutoff in the field theory 
analysis. Our choice of U and V makes the numerical analy- 
sis somewhat easier and shows that the Ising quantum phase 
transition is a robust feature of the lattice model. 



A. Order Parameters 

The DMRG calculations of the order parameters 0, Q 
were performed in chains with open boundary conditions 
(OBC). We used up to L = 1024 lattice sites and kept up to 
to = 1024 density-matrix eigenstates in the truncation of the 
superblock Hamiltonian. The results are summarized in FigQ 
At zero dimerization, the system is a CDW-insulator for any 
value V > U/2 12111 . When we turn on the dimerization, 
the BOW order parameter (tobow) grows linearly with 6 in 
agreement with the semi-classical analysis d!61 i and the per- 
turbative result ( 1241 . While (tobow) is enhanced, the charge- 
density wave parameter (tocdw) is reduced until it decays 
rapidly at the quantum critical point 



rOBC 



1.28. 



(26) 



Beyond this point (tocdw) vanishes, whereas (tobow) devi- 
ates non-trivially from the line (tobow) = bo8 with bo(U jt — 
4,V/t — 3) w 0.25. For small values of S we can see that 
(tocdw) = a o — a i<5 2 for some constants a,Q and oi. This 
is in agreement with the predictions d23i and d!6i . The onset 
of (tocdw (8)) as a function of S close to the critical point 
is strongly reminiscent of the magnetization in the classical 
two-dimensional Ising model. We therefore attempt a fit of 
the data with a power-law onset (tocdw(^)) ~ c o|^ — ^cC 3 in 
the vicinity of the critical point and find (3 = 0.126 w 1/8. 
The logarithmic plot in the inset of figure ^ shows the good 
agreement of this fit with our data. 
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FIG. 1: Bond-order wave ((meow)) and charge-density 
wave ((meow)) order parameters of the extended Peierls-Hubbard 
model with U/t = 4, V/t = 3, and varying dimerization 8. There is 
a sharp transition at 5® BC — 1.28. For small dimerization 5 clearly 
(mBOw(<5)) oc S (full line). Inset (a): The CDW order parameter 
diminishes quadratically. The dash-dotted line is a quadratic fit of 
the form i 1 61 . Inset (b): log-log plot of the CDW order parameter 
and a fit with the power-law onset (iticdw) ~ \S — 5 C | 0,126 (dashed 
line). 



This confirms our suggestion that the transition belongs to 
the Ising universality class. A finite-size scaling analysis of 
spin and charge excitation gaps in the following subsection 
further corroborates this conclusion. 



B. Excitation Gaps 

We define the spin and one-particle gaps 

A s = E (N,1) - E (N,0) , 
A c = E (N + 1, 1/2) + E (N- 1,1/2) 
-2E (N,0), 

and the gaps to the first and second excited state 



Ai = E 1 (N,0)-E o (N,0), 
A 2 = E 2 (N,0)-E (N,0). 



(27) 
(28) 



(29) 
(30) 



In these definitions E (N, S z ) is the ground state energy in 
the subspace with a given number N of electrons and a given 
spin S z . Likewise, E\{N, S z ) and E%(N, S z ) are the energies 
of the first and second excited state, respectively. 

In contrast to the previous section, we do not employ open 
boundaries to calculate the excitation gaps since we find that 
localized bound states occur at the system boundaries. Since 
we are not interested in the energy of such surface states, we 
use periodic boundaries (PBC). We studied periodic chains 
with an even number of lattice sites and chain lengths up to 
L = 128 while keeping as much as m = 3072 density-matrix 
eigenstates. 

Figure|2]shows a plot of the gaps as a function of the dimer- 
ization. The gaps Ai and A 2 are strongly size dependent. 



Therefore, we apply finite-size scaling analysis to extrapo- 
late the gaps to the thermodynamic limit. Below the critical 
dimerization the gap Ai extrapolates to values very close to 
zero. This means that the ground state is twofold degenerate 
in the CDW phase. Above the critical dimerization the gap 
Ai opens linearly and the ground state is no longer degen- 
erate and displays no long-range CDW order. At S = the 
extrapolated gap to the second excited state A 2 is very close 
to the value of the spin gap A s which we expect to be equal 
in a CDW-insulator. They stay close also for small dimer- 
izations which indicates that the CDW phase of the extended 
Hubbard model is not strongly perturbed by a small dimeriza- 
tion. Tuning 6 to larger values the spin gap A s is not much 
affected in contrast to A 2 which is now linearly reduced with 
growing dimerization. Above the transition A 2 increases with 
the dimerization. Figure|2]suggests that A 2 is at most slightly 
larger than Ai in the thermodynamic limit or possibly degen- 
erate. We find that the one-particle gap A c > A s for any 
dimerization. 
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FIG. 2: Dimerization dependent gaps Ai, A2, A a , and A c . In the 
interval S < 8 C the ground state is degenerate and Ai = 0. The 
gap A2 is reduced linearly (full line) as we approach S c . Above the 
transition the gap to the first excited state, Ai, opens linearly (dashed 
line). The spin gap A s is non-zero for all S and is equal to A2 for 
very small dimerization. We always have A c (<5) > A s (<5) for the 
one-particle gap A c . 

In order to determine the critical point for the periodic 
chains we calculated the dimerization dependence of A 2 for 
many different system sizes. The results are shown in the in- 
set of figure |3 We determine the minima min[A 2 (5)] and 
their positions S m [ n (L) by fitting second order polynomials 
to the curves A 2 (J) for various system lengths L. We then 
extrapolate these quantities to the thermodynamic limit. This 
is shown in figure [3] where we observe that A 2 (J m ; n ) — > 
within the precision of our extrapolation. The value of the 
critical dimerization determined from S m i n (l/L = 0) reads 



?PBC 



= 1.29, 



(31) 



in good agreement with the result ( I26> previously obtained 
with open boundaries. 
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FIG. 3: Extrapolation of the minima min[A2 (<5)] and their positions 
<5min to the limit 1/L — > for periodic boundary conditions. Within 
the numerical precision min[A2(l/Z/ —* 0)] = and BC = 1.29. 
Inset: A2(<5) for various L. 

We have seen that the onset of (tocdw) is compatible with 
an Ising-type phase transition. Now, we can go further to show 
that the excitation of the system that becomes critical at S c also 
suggests this interpretation. As 8 approaches 8 C we expect that 
the gap to the lowest excitation vanishes like 

A ± ~A ± \S-6 C \ ZV , (32) 

below (— ) and above (+) the critical point. The non-universal 
constant is a typical energy scale of the system and zv 
is a universal critical exponent. In figure [2] we show a linear 
fit of the extrapolated gap to the lowest excited state above 
(A + = Ai) and below (A - = A2) the transition point. Both 
fits are a good description of our data since we can derive a 
critical dimerization 5 C « 1.25 which is consistent with (126 \ 
and OH . We can now infer that zv = 1. In order to fix z we 
note that the characteristic length scale £(8) of the critical fluc- 
tuations diverges at the critical point such that £ ^ |<5 — ^ 



holds. The length scale £(8) can be estimated by consider- 
ing the critical dimerization 8 C (L) as a function of the system 
length L. By inverting this relation we obtain a critical system 
size L c (8) which is an estimate of the length scale £(<$) of the 
critical fluctuations. From this we find that v = 0.98 w 1. 
Since the characteristic energy scale A* vanishes linearly we 
conclude that the dynamical critical exponent 2=1. Since 
both (3 w 1/8 and v w 1 are independent universal exponents 
we may conclude that the observed quantum phase transition 
belongs to the Ising universality class. 



IV. CONCLUSIONS 

We have shown that there is a quantum phase transition 
from a mixed CDW-BOW to a BOW phase in the half-filled 
extended Peierls-Hubbard model. Field theory arguments 
suggest that this phase transition belongs to the universality 
class of the two-dimensional Ising model. A DMRG study of 
the extended Peierls-Hubbard model for parameters U/t = 4, 
V/t = 3 reveals that there is a transition at a critical value 
8 C ~ 1.3, where the CDW order parameter (mcDw) is found 
to vanish. A detailed analysis of the order parameters and ex- 
citation gaps in the vicinity of the transition confirms that the 
transition falls into the Ising universality class. This is a robust 
property of the lattice model away from the weak-coupling 
limit. 
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